WEAK ENTROPY SOLUTIONS OF NONLINEAR 
REACTION-HYPERBOLIC SYSTEMS FOR AXONAL TRANSPORT 



HAO YAN AND WEN- AN YONG 

Abstract. This paper is concerned with a class of nonlinear reaction- hyperbolic sys- 
, terns as models for axonal transport in neuroscience. We show the global existence of 

entropy-satisfying BV-solutions to the initial-value problems by using hyperbolic-type 
methods. Moreover, we rigorously justify the limit as the biochemical processes are 
much faster than the transport ones. 

Keywords: axonal transport, relaxation limit, difference scheme, BV-estimates, en- 
tropy. 
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1. Introduction 



The axonal transport is important for the maintenance and functions of nerve cells. 
These cells are also called neurons. A neuron consists of three parts mainly: cell body, 
dendrites and a single axon. The axon is a long and thin pipe whose length can exceed 
10,000 times its diameter. It is this axon that distinguishes neurons from other cells. The 
axon is responsible for signal transmission in the nervous system. Its cytoplasm does not 
contain rough endoplasmic reticulum and therefore its proteins can only be transported 
from the cell body, where all proteins are synthesized. 

The transport proceeds as follows. Proteins are stored in vesicles as cargos. The vesi- 
cles are attached to kinesin (anterograde motors) or dynein (retrograde motors) proteins. 
These motor proteins drive the vesicles to walk along the cytoskeletal microtubules as 
track. Here the kinesin proteins move the vesicles from the cell body to synapse (antero- 
grade transport), while the dynein proteins move the vesicles in the opposite direction 
(retrograde transport). During the transport, many biochemical processes are possible. 
For example, the cargos can leave its track, can switch its motor proteins from kinesin 
to dynein or vice verse, and can move back onto the track. Thus, we can divide the 
cargos into a number of subpopulations, such as free vesicles, vesicle-kinesin compounds 
off track, moving vesicle-dynein compounds on track, etc. 

l 
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As the axon is long and thin, it is reasonable to assume the transport only along the 
longitudinal direction of the axon. Denote by x > the distance down the axon from the 
cell body which is located at x — 0. Let Ui = Ui(x,t) be the concentration at space-time 
(x, t) of the i-th subpopulations. According to Reed and Blum [6] , the mathematical 
model for axonal transport is partial differential equations of the form: 

d t Ui + \id x Ui = F i (u 1 ,u 2 , ■■■ ,u r ), i = 1, 2, • • • , r. 

Here the term \id x Ui accounts for the transport of the z-th subpopulation with constant 
velocity Aj, and F t (ui, u 2 , ■ ■ ■ , u r ) describes the biochemical processes of the constituents. 
It is well recognized that the biochemical processes are much faster than the transport in 
biosystems. 

In [21 IH El [7], the authors studied the linear case, where Fi(ui,U2, ■ ■ ■ ,u r ) is linear 
with respect to the u/s, in order to explain the approximate traveling waves observed 
in experiments. Especially, in [5] Friedman and Hu used parabolic-type estimates to 
analyse the diffusive limit of the linear systems. However, it seems uneasy to deal with 
the nonlinear problems with the parabolic-type techniques. On the other hand, in pQ 
Carr showed the existence of global classical solutions to a class of nonlinear models with 
source terms of the form 

J7i, i = i 

F i (u 1 ,u 2 , ■ ■ ■ ,Ur) = <fi - fo-i, 1 < i < r , 
[-fr-i, i = r 

where /j = fi(v,i, u i+ i) is a continuously differentiate function of two variables. We notice 
that in applications fi is a polynomial of itj and Uj+i. In addition, the model in [2] is an 
example where /, depends only on and u i+ i. 

In this paper, we consider the same nonlinear systems, as in [lj, but with a small 
parameter e > 0: 

(dt + XidJinfat) =-f 1 {ui,u 2 ), 



(d t + Xid x )ui(x,t) = --fi^Ui-^Ui) + \fi{u h u i+1 ), 



(d t + X r d x )u r (x,t) = -~f r _i(u r -i,u r ). 

Here the small parameter e characterizes the fact that the biochemical processes are much 
faster than the transport. We assume, throughout this paper, that each fi is strictly 
decreasing with respect to the first argument and strictly increasing with respect to the 
second. In addition, we assume that /j(0, 0) = and, for fixed v, there exists w such that 
vfi(v, w) > 0. These assumptions are consistent with those used in PQ [21 111 El El E] • 

We will regard (11. ip as a hyperbolic relaxation system [10] and use the corresponding 
techniques to study it. In particular, we will use a difference scheme to show the global 
existence of entropy-satisfying BV-solutions (u\,u 2 , ■ ■ ■ ,uf) to the initial- value problems 
of (11.11) and investigate the limit as e goes to zero. For r < 3 and for linear problems, 
we prove that the limit (u1,u 2 , ■ ■ ■ ,u®) is also an entropy-satisfying BV-solution to the 
so-called equilibrium or reduced system of (II. ip : 

d t (ui + u 2 -\ h u r ) + d x (XiUi + X 2 u 2 H h X r u r ) = 0, 

fl(Ul,U 2 ) = /2O2, = • • • = /r-lOr-l, «r) = 0. 
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Although the BV-frsanework is quite standard for nonlinear hyperbolic problems (see, 
e.g., [Hill E]), some innovatory ideas seem necessary to carry our the details. In particular, 
we use the Brouwer fixed-point theorem, the i?^/-estimate and the special structure of the 
source terms in (11 .ip to derive the existence, boundedness and time-Lipschitz continuity 
of the difference solutions. It seems not so easy to obtain the time-Lipschitz continuity! 
Moreover, for an arbitrarily given convex entropy function for the reduced system (11.21) . 
we construct a dissipative entropy function for the original system (II. ip . In equilibrium, 
the constructed entropy function reduces to the given one. 

The paper is organized as follows. Section 2 is devoted to the analysis of a difference 
scheme for the system (11.11) . The time-Lipschitz continuity of the difference solutions is 
derived in Section 3. In Section 4, we discuss entropy functions for the two systems (11.11) 
and (11.21) . The main results are shown in Section 5. 

2. Difference Solutions 

We begin with construction of approximation solutions to the reaction-hyperbolic sys- 
tem (II .ID by using difference methods. For simplicity, we set 

U= (ui,u 2 , ■ ■ ■ ,u r ) T , A = diag(Ai, A 2 , ■ ■ • , A r ), 

Q(U)= (/!,••• ,-fi-x + fi,--- ,-/^i) T , 

where the superscript T denotes the transpose of vectors or matrices. Then the reaction- 
hyperbolic system (II. ip can be rewritten as 

(2.1) U t + kU x = -Q(U). 

About this system, we make the following assumptions mentioned in the introduction: 

(1) . each fi is continuously differentiable, strictly decreasing with respect to the first 

argument, and strictly increasing with respect to the second; 

(2) . £(0,0) =0. 



These assumptions are consistent with those used in [IJ EJ HI El El E] • 

Our difference approximation to (12. ip is the following semi-implicit upwind scheme 

jjn+l _ Tjn Tjn _ Tjn Tjn _ Tjn , 



At Ax ~ Ax e 



Here Ax and At denote the increments respectively in x and t; denotes the approxima- 
tion of U (x, t) over the grid block [xj,Xj + x) x [t n , t n+ i) with Xj = j Ax and t n = nAt; j = 
0, ±1, ±2, • • • and n = 0, 1, 2, • • • ; Q] = Q{U?); A + = diag{Xf, A+, • • • , A+) and A" = 
diag(Xi, A 2 , • ■ ■ , A~) with 

(2.3) K = ^^>0 and Ar = < 0. 
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For n = 0, we take 

(2.4) U°=l 

U if j>-k 

where Uq(x) := (uio(x), U2o(x), u r o(x)) T is a bounded measurable function of x G R. 
In addition, the in (j2.4p can be replaced with any constant vector. 

Throughout this paper, we assume that the grid sizes satisfy the CFL-condition 

2.5 — max Ai < 1. 

( J Ax i 1 1 ~ 

To analyse the above scheme, we start with the following elementary fact. 

Lemma 2.1. Let A = (cifcz)nxn be a matrix satisfying X)fc=i a ki = 1 f or I = 1)2, ■ •• , n 
and aw — for k ^ I. Then A is invertible, the spectral radius of A^ 1 is not bigger than 
1, and the 1-norm of A' 1 is 1. 

Proof. Let A be an eigenvalue of A. By the Gershgorin circle theorem there is an integer 
/ such that 

|A - a u \ < 2j \a k i\ = an - 1, 

k^l 

for an — 1 — J2k=a a ki an d a ki < 0(fc 7^ /). From this we easily deduce that |A| > 1. Thus 
A is invertible and the spectral radius of A -1 is not bigger than 1. 
Furthermore, it follows from ^2=1 akl = ^ ^hat 

(1,1,...,1)A=(1,1,...,1) 

and thereby 

(1,1,... ,1)A~ 1 = (1,1,--- ,!)■ 
This shows that the sum of each column of A" 1 is also 1. On the other hand, we set 
D = diag(an, a 2 2, • ■ ■ , a nn ) and F = D — A. It is obvious that the elements of FD~ l are 
nonnegative, 

A = D-F = (I-FD-^D 

and thereby 

A' 1 = D~ X {I - FD' 1 )- 1 . 

Moreover, with the Gershgorin circle theorem it is not difficult to see that the spectral 
radius of FD~ X is less than 1. Thus we can write 

A' 1 = D-\I + FD- 1 + (FIT 1 ) 2 + ••■), 

which shows that the elements of A~ x are nonnegative. Hence the 1-norm of A" 1 is 1 for 
the sum of each column of A' 1 is 1. This completes the proof. □ 

Set V = (v 1, v 2 , ■ • • , v r ) T , W = (wi, w 2 , ■ • ■ , w r ) T , 

(2.6) At = A(V,W)= [ ^(9v i + (l-e)w i ,ev i+1 + (l-9)w i+1 )d9, 

Jo OUi 

(2.7) Bi = Bi(V,W)= f -^(9v i + (l-9)w i ,ev i+1 + (l-e)w i+1 )de. 

Jo ° U i+l 
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Assumption (1) implies that Aj < and Bi > 0. By the mean-value theorem we have 
(2.8) Q(V) - Q{W) = Q(V, W)(V - W). 

From the definition of Q, we have 

/ A x Bi 



Q{V, W) 



\ 



-At -B x + A 2 B 2 

-A 2 -B 2 + A 3 



B, 



V 



-A._l -Br-! J 



It is easy to verify that / — ^Q(V, W) satisfies the conditions in Lemma [2.11 Therefore 
we have 



Corollary 2.2. / - ^Q(V, W) is invertible, the spectral radius of (I - ^Q(V, W))' 1 
not bigger than 1, and the 1-norm of (I — — Q(V, W))' 1 is 1 for any V, W E M r . 



IS 



With these preparations, we study the difference scheme (12.21) . In what follows, we will 
denote by | ■ | the 1-norm of vectors and matrices. 

Lemma 2.3. Given 11™, U™ +1 can be uniquely determined by solving the nonlinear alge- 
braic equation (\2.2S) . 

Proof. From fi(0, 0) = we see that Q(0) = and 

Thus, we deduce from the difference scheme (12.21) that 
At 

:;(/■■'"-' nn/''" n_ = — r r / " +(/- — v r - 

Ax 



rn+l 



(I 



-Q(f/ 7 n+1 ,0))f^ 



Ax Ax 



With this equation, we construct a mapping F : V — > W as follows 



(/ - *Q(y, Q))W = ^ + uu + (/ - 

Due to Corollary 12. 2\ we can write 

W = (I- ^QivM-'^uu + (I - ^a- 



At 
Ax 



A-)f/jVi- 



Ax ' 3 v Ax ' 3+11 



and thereby 
\W 



< |(/ - ^Q(V, O))-!!^^ + (/ - £ A + + ^A~)U- + (-^A-)^xl 



Ax v Ax Ax 



'Ax 3 Ax Ax ' 3 v Ax 

Therefore, W is bounded for fixed j and n. By the Brouwer fixed-point theorem, the 
mapping F has fixed points for F is continuous. This shows the existence of . 

To see the uniqueness, we assume that both U\ and C/2 satisfy the difference scheme 
02. 2p . With the mean- value theorem, we see that 



Ut-U 2 



At 



{Q{Ui)-Q{U 2 )) 



At 



■Q{U 1 ,U 2 )(U 1 -U 2 ) 
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and thereby 

(I-^Q(U 1 ,U 2 ))(U 1 -U 2 ) = 0. 

It follows from Corollary 12.21 that U\ = Ui- This shows the uniqueness and hence the 
proof is complete. □ 

Remark 2.1. We could use the contraction mapping principle to show the existence and 
uniqueness of . However, it requires that At is much smaller than e. On the other 
hand, this lemma does not tell the uniform boundedness of the difference solution 

Next we establish the L 1 -stability of the difference scheme (12.21) . 

Lemma 2.4. Let U™ and VJ 1 be two solutions of the difference scheme (\2.2l) with initial 
data Uj and V®, respectively. Then it holds that 

+oo +oo 

£ w- v r\< £ W-W 

j = — OO j = — OD 

for all n > 0. 

Proof. It follows from the difference scheme (I2.2p that 

Tjn+i_ V n + i = ^(g(E^+ 1 )-g(^ 1 )) + ^A+(t^_ 1 -^ r L 1 ) 

Ax Ax J J Ax J J 

Thus we deduce from the mean- value theorem and Corollary 12.21 that 
At At 

(/ - —Q(u? +1 , v; i+l ))(u] i+l - v; +1 ) = ^ A+ (^ n -i - y£0 

jjn+l _ yn+1 = (J _ ^Lq{JJ^\ V^ 1 )) ~ X A + (C/J l _ 1 - VJLj) 

+V - |^A+ + ^A-)(U- - V?) + (- £ A")(^ +1 - V? +1 )], 
\U? +l -V 3 n+1 \ < \(I-—Q(U? +1 ,VJ l+1 ))- 1 \\^A+{U?_ 1 -VJ l _ 1 ) 



< 



e x 3 ' 3 " "Ax 
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and thereby 

+00 +00 . 

E wr i -vr\< e itw a+ (^-^)i 



+°° +°° A + 

+ E i( / -^ A+ + ^)( f/ r-^)i+ E k-^-)W + i-*k 



hoo 



IS 



j=-oo 

+I(-^a-)(c/;-^)|]. 

In the last step we use the induction assumption that X]j=^oo I^T — ^Tl < 00 > which i 
true if X^woo — < 00 • On ^ ne °ther hand, from the definition of (I2.3P and 
the CFL-condition (12.51) we deduce that 

\^(u« - v;)\ + \(i - |^a+ + £a-)V? - yp)\ + K-^A-Xi? - 
= E ^ W - «5I + E(i - + - «SI - E |^K- - «3I 

i=l i=l i=l 

= [[/!> _ y»| 

Hence we see that 

+00 +00 +00 

E i^ +1 -^ +1 i< E w-v?\<...< e w-^ri- 

J= — OO j = — OD j = — CO 

This completes the proof. □ 

By taking VJ 1 = £/]Li in Lemma [2. 4l we get the following corollary on BV-estimates of 
the difference solutions. 

Corollary 2.5. Let U™ be a solution to the difference scheme (\2.2l) with initial data If®. 
Then the BV -estimate 

+00 +00 

E W-uu\< E M- U U\ 

j = — OD j = ~OC 

holds for all n > 0. 



Having Corollary I2.5I we show the uniform boundedness of the difference solutions, 
which is not covered in Lemma I2.3I 

Lemma 2.6. 

+00 

sup|t/?|< E W-^-il 
3 

3=-oo 

holds for all n > 0. 
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Proof. Thanks to the initial data ( 12. 4p . we have U™ = if jo is large enough. Thus, for 
any j it follows that 



.70 



\u?\ = \u?-u%\< \Uk-uiU\< \ u k- u °k-i\- 

k=j+l 

This completes the proof. 



k=— oo 



□ 



Remark 2.2. This lemma requires that the initial data are of bounded variation. Other- 
wise, we could not obtain the uniform boundedness of the difference solutions. 



3. Time-Lipschitz Continuity 

In this section we show the time-Lipschitz continuity of the difference solutions, which 
seems not so easy for the present problem. 

To begin with, we set G(U) = (A, / 2 , • • • , f r -i) T G M r_1 and G] = G(U] 1 ). Note that 

(3.1) Q(U) = KG(U) } G(U) = K'Q(U), 

where K is the constant r x (r — lVmatrix 



/ 1 



(3.2) 



K 



\ 



-1 1 



V 



-1 / 



and K' is the constant (r — 1) x r-matrix such that K'K = I r -i- Moreover, it follows 
from the mean-value theorem and the definitions in (12.61) and (12. 7\i that 



G n+1 _ G n = M (Jjn^ 



v?) 



with MiU^ £/ n+1 ) the following (r — 1) x r matrix 



M{U?,U? +1 ) 



( A l B l 

A 2 B 2 

A 3 B 3 



\ 



\ 



A--1 B r _i J 



Set M = M{U?, U™ +1 ) = M(Uf, U™ +1 )K e R(r-V*<r-i) t By ca i cu i at ion, we have 



M 



( Ai - Si Bi 

— A 2 A 2 — B 2 B 2 

—A 3 A 3 — B 3 B 3 



\ 



-A r _i A 



r-l 



\ 



B r -i J 



Thanks to the strict monotonicity assumption (1) on the /j's, it is an elementary fact that 
I _ £±m is invertible. 
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We assume that there exists a ([/", £/™ +1 )-independent norm | • |* on R/*" 1 and a positive 
constant A such that 

(3.3) M)- 1 U<(1 + — A)" 1 . 

e e 

Although such a norm has not been found for general cases with r > 3, this assumption 
is indeed true for r < 3 and for the case where Q(U) is linear with respect to U. In fact, 
we have 

Proposition 3.1. For r < 3, there exists a positive constant A such that 

|(/-— Af)- 1 |oo<(l + — A)" 1 , 
e e 

Proof. When r = 2, M is the negative number (Ax — B{). Thus we have 

|(/ - — M)-^ = (1 + — (Si - AO)" 1 < (1 + —A)- 1 
e e e 

with A = miny i H/{-Bi(^ // ; W) — Ai(V,W)} > 0. Here the boundedness of the different 
solution established in Lemma 12.61 has been used. 
For r = 3, M is the 2x2 matrix 



M 

Then we have 



-A 2 A 2 - B 2 



and 



I _At M _ l + ^B.-Ax) -f Bl 
1 e M ~{ fA 2 l + f(B 2 -A 2 ) 

1 (i + ^B 2 -A 2 ) M Bl 

l + A ^ + B ^L\ -fA 2 i + M( Bl - Al 

with A = Bx - A x + B 2 - A 2 and B = B X B 2 - A X B 2 + A X A 2 . Recall that A, L < and 
5,; > for i = 1,2. Set 

a = max^! + B 2 - A 2 , B 1 -A 1 - A 2 }. 
It is obvious that < a < A and B > 0. Thus we deduce that 

with 

A = min{A - a,B/a} > 0. 
This completes the proof. □ 

For r > 3, M is no longer strictly diagonally dominant and the inequality in Proposition 
13.11 does not hold anymore. However, we have 

Proposition 3.2. If Q(U) is linear with respect to U, then the above assumption (\3.3\) 
holds. 
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Proof. Since Q(U) is linear with respect to U, M is a constant matrix. Set a\ — 1, ctj + i = 
— ttj^^-, and .D = diag(ai, a 2 , ■ ■■ , a r _i). It is obvious that the diagonal matrix D is 

positive definite and DM is symmetric. Then D^MD~2 is also symmetric. Thus, there 
exist an orthogonal matrix T and a diagonal negative-definite matrix N such that 

D^MD~^ = T~ l NT, 

and thereby 

(TD*)(I - ^M)- 1 = (I- ^N)- l {TD^). 
Set P = TD^. We define a norm on R r_1 as 

\z\ P = \m 

for £ € R r_1 . Thus the corresponding matrix norm is 

|(/-— MY\ := sup \(I-—M)- 1 a P 
e \S\p<i e 

= sup \P(I-—M)-^\ 

I€Ip<i e 
At 

= sup |(J AT)- X P£| 

I€Ip<i e 

< sup |(/_^iV)-i||P£| 
KIj»<i e 

< |(J - ^JV)" 1 ! = (1 + ^A)" 1 , 

where A is the smallest eigenvalue of — N. □ 

Now we turn to estimate Q n = Q{U^). 

Lemma 3.3. Assume hS.ty holds. Then there is a positive constant C such that 

\\Q n \\ L ,<C((l + \^)- n \\Q°\\ L i+e). 

for all n. 

Proof. From the difference scheme it follows that 

G n+l_ G n = M(U^U] +1 )(U™ +1 -UJ) 



and thereby 



At 



At 



+A-(c/; +1 - c/;)] 



Ax' 
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Thus we deduce from the inequality (I3.3P that 



A*,,^ At 



1 + A— < |(/ - — M)G n +\ 

At 



e 

\G] - M(U?, U? +1 )[A + (U? - U?_ x ) + A"(E^ +1 - Uf)} 



Ax 



At 



< \G]U + \M(u?,u? +1 )[A+(u? - uf_ x ) + a-(c/; +1 - c/jv Ai . 

Since | • |* is equivalent to the 1-norm | • | on R r_1 and since the difference solution is 
bounded, there is a positive constant C such that 

(i + \^)\GTX < \G]U + c£(\u? - U^\ + \U» +1 - U»\). 

Consequently, we use Corollary 12.51 to obtain 

(1 + A— )^2\G]+%Ax < \ G j\* Ax + CAt ^2\ U j ~ u j-i\ 



< J2\G]UAx + CAt. 



From the last inequality, it is easy to verify that 

\G]\,Ax < (1 + A—)"" \G%Ax + Ce< C((l + A— )~ n ^ Ax + e) 

3 3 3 



Thus we get 



|G n || L i < Cj2\ G ]\* Ax < C((l + A— )- n ^|G°|Ax + e). 

3 3 



Finally, from (13.1 ft it follows that 

\\Q n \W <C\\G n \\ L , and ||G°|| L i <C||Q°|| L i. 



Hence 



\\CF\W<C{{\ + \— )- n \\Q°\W + e) 
e 



and the proof is complete. □ 

Now we can easily show the time-Lipschitz continuity of U n . 
Lemma 3.4. Assume A3. 3\) holds. Then 

- w\\ L i < c(— (i + a— )- (n+1) ||g || L i + At) 

for n = 0, 1, 2, • • • . 

Proof. From the difference scheme ( 12. 2 p we have 
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and therefore 

Wr 1 - u-\ < ^\Q] +1 \ + c£{\u? - u?_ x \ + \v» +l - u»\). 

Thus, from Corollary 12.51 and Lemma 13.31 we deduce that 

\m +1 - U?\\» < — HQ7+ X \W + CAt < C(— (1 + A— )- (n+1) ||g°|| Ll + At). 
6 e e 

This completes the proof. □ 



4. Entropies 

In this section, we discuss entropy functions for the reaction-hyperbolic system f 1 1.1 1) 
and its equilibrium system (11.21) . For this purpose, we make the following additional as- 
sumption 



(3). For any fixed v, there exists w such that vfi(v, w) > for i = 1, 2, • • • , r — 1. 



This assumption, together with those made in the previous section, ensures that there 
exists a unique and globally-defined function hi of one variable such that 

(4.1) fi(ui,u i+ i) = iff u i+1 = hi(ui) 

for i — 1,2, ■ • ■ ,r — 1. Obviously, hi is strictly increasing and hi(0) = 0. By the implicit 
function theorem, hi is continuously different iable. Thus the equilibrium system (11.21) can 
be written as 

d t (ui + u 2 -\ h u r ) + dJXxUi + \ 2 u 2 H V \ r u r ) = 0, 

(4.2) 

«2 = M 3 = h 2 (u 2 ), • • • , Mr = /ir-l(Wr-l)- 

Set 

v — ui + hi{u\) + h 2 o /i 1 (ni) + • • ■ + h r _i o /i r _ 2 o • ■ ■ o ^i(«x)- 

Since the right-hand side is strictly increasing with respect to u±, U\ can be expressed as 
a function of v, say Ui = ui(v). By the inverse function theorem, Ui(v) is continuously 
differentiable. Set 

u i+1 (v) = hi(ui(v)) 

for i — 1, 2, ■ ■ • , r — 1. Consequently, each Ui(v) is strictly increasing, continuously differ- 
entiable and Ui(0) = 0. Set 

v = U\{v) + u 2 {y ) + • ■ ■ + u r {v ) and h{y ) = XiUi(v) + X 2 u 2 (v) + ■ ■ ■ + \ r u T {v). 

Then the equilibrium system (14.21) can be rewritten as 

(4.3) d t v + d x h{v) = 0. 

Recall that any convex function is a convex entropy function for scalar conservation laws 
like (USD (see, e.g., 0). 

Next, we turn to discuss the entropy functions for system (14.31) and the reaction- 
hyperbolic system (II. ip . 
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Lemma 4.1. Given a strictly convex smooth function fj(v), there is a dissipative entropy 
function r)(U), in the sense of [IT] , for the reaction-hyperbolic system fli.il) such that 

fj{v) =T](u 1 (v),U 2 (v),--- ,U r (v)) =T](U(v)). 

Proof. We inductively define 



T} r (u) = / fj(u r (w))dw, 
JO 

?7i_i0) = / T]' i (h i -i(w))dw 
Jo 

for i = r, r — 1, • • • ,2. Then we have 

(4.4) Vi-M-i) = iiih-iiui-t)) and r)' r (u) = v'(u;\u)). 

Since 77(f) is strictly convex and u~ x and hi are all strictly increasing, it is clear that 
771,772, • •• ,77^. are strictly increasing. Therefore, 

r 

(4.5) 77(C/):=^77,( Ui )+77(0) 

is a strictly convex function of U. Recall the matrix K defined in (13. 2p . We see that 

and, moreover, from (13. IJ that 
Q(U) = KG(U) 

(A a\ ISA- //l( M l' M 2) f2(u 2 ,U 3 ) fr-l(u r -l,U r ) T 

(4.6) = ATdiag(— — , — — , • • • , —)K -quiU) 

Vl -V2 % - % Vr-l - Vr 

= S(U) Vu (U). 

Thanks to the relations in (14. 4p . we deduce from (14. ip . the convexity of 77^ and the mono- 
tonicity assumption (1) that 

fi{Ui,U i+1 ) fi(Ui,U i+1 ) 



Jo fiu i+1 (ui, hi{ui) + a{u i+1 - hi{ui)))da 



< 0. 



Jo Vi+i( h i( u i) + a ( u i+i - hi{ui)))da 
Thus, S(U) is a symmetric and non-positive definite matrix. Its null space is obviously 
that of K T , which is independent of U. Consequently, rj{U) is a dissipative entropy 
function, in the sense of [UJ, for the reaction-hyperbolic system (11.11) . 
Furthermore, we deduce from the relations in (14. 4p that 

v[{M v )) = v' 2 ( h i( u i( v ))) 

= ^{h 2 o hxiux)) 



r)' r {h r -i o h r _ 2 o • • • o hi(ui)) 
r)' r {u r {v)) 
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Similarly, we have 
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T]' 2 (u 2 (v)) = ■■■ = T]' r (u r (v)) = fj'(v). 



Therefore, we deduce from U\{v ) + u 2 (v ) + •■■ + u r (v) = v that 

(v(u(v)))' = ($>(^)))' 

8=1 

= ViM v )) u 'i( v ) + V2{u2(v))u' 2 (v) H h rf r (u r (v))u' r (v) 

= V'{v). 

In view of rji(0) = = Mj(0), we see from (14.51) that 1](U(0)) = rj(0) and hence 

fj( v ) = T]{U(v)). 

This completes the proof. □ 
We conclude this section with a discrete entropy inequality for the difference solutions. 

Lemma 4.2. Let U r f be a solution to the difference scheme f lff.ffl) . Then, for any smooth 
convex function r](U) = Y7i=i Vi( u i)> there exists a Lipschitz continuous function ^ of two 
variables such that for all jeZ and n > 0, the following cell entropy inequalities hold: 

rn+l\ ^ ^/TTn\ i ,t, / j jn jjn \ \\i/TT n TT n \\ i (TT nJ r^- \ f^tl T /""■+! '\ 



j 



Moreover, the Lipschitz continuous function satisfies the consistency relation 

r 

*(U,U) = £ 
i=i 

Proof. For any smooth convex function rji and any two real numbers a, b G R, it is standard 
that 

r)i(b) - r)i(a) < 77 ■(&)(& - a). 
Thus, for the given convex function r](U) = Yli=i Vi{ u i) an d for any V, W e W we have 

rj(V) - r){W) < r)u{V){V - W). 
Thus, it follows from the original difference scheme that 

V(U? +1 ) < V(U? - |^A + (f/« - U?_ x ) - ^-(C/JV! - U^)) + Vu (Uf)^Q^\ 
since 
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On the other hand, we deduce from (12.31) . (I2.5P and the convexity of r\i that 



Ax 7-1 v Ax Ax ' J Ax 
I^tf^-i + ^Ax + Ax ^ - * Ax"<^ ' 



i=l 

r At 

i=i 



Define 



i=l 

This \I> is obviously Lipschitz continuous and satisfies the consistency relation. Moreover, 
the above inequalities lead directly to 

This completes the proof. □ 



5. Main Results 

In this section we prove the main results of this paper. To this end, we define 

A/„ l\ /„.A/ ,\ A/ ,\ „.A/ ,\\T . /„ ,n „.n „,n \T 



U (X,t) = (ttj (x,t),tt 2 (x,t), ■ ■ • ,M r (x,t)) := « >i ,'U2j 



for (x, t) G [jAx, (j + l)Ax) x [nAt, (n + l)At). With this definition, it simply follows 
from Corollary 12.51 Lemma 12.61 and Lemma 13.41 that 

Lemma 5.1. The piecewise constant function U A (x,t) satisfies the following estimates 

(5.1) \U A (x,t)\< \ U !~ U j-i\ forall(x,t), 

j=-oo 

+ 0O 

(5.2) TV(U A (-,t))< \Uj-Uj-il 

j=-oo 

(5.3) \\U A (;t) - f/ A (-,ti)|| Ll < C(- e (l + A^)-"^||g || Ll + l)(|t - h\ + At) 
for all t, t\ > 0. 
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Proof. We only need to show the last inequality. Let k, k\ be two integers such that 
t G [kAt, (k + l)Ai),ti G [kiAt, (k\ + l)At). Without loss of generality, we assume 
fci < k. Then we deduce from the definition of U A and Lemma [3.41 that 

||^(., £ )-^A ( . )fi) | Ui = J2\U*-Uf\Ax 

j 

n=ki j 

< Vc(-(i + a— )-(" +1 )||g°|| il + i)At 



n=k\ 

< C(-(l + A— r (fcl+1) ||Q°|Ui + l){k - h)At 

e e 

1 At min-ft.t-i } _ 

< c(-(l + A— *niQ°|| L i + 1) |t - til + At). 

e e 

This completes the proof. □ 



Having the estimates in Lemma 15.11 and the discrete entropy inequality in Lemma 14.21 
we follow the standard argument in [3j to obtain (e is fixed) 

Theorem 5.2. Suppose the initial data Uq(x) = (u w (x), U2o(x), ■ ■ ■ , u r o(x)) have bounded 
variations, the grid sizes At and Ax satisfy the CFL-condition (\2.5h . the f^s satisfy the 
assumptions (l)-(3), and the inequality hS.^I holds. Then, as the grid sizes At, Ax tend to 
zero, there is a subsequence of 'the function family U A (x,t) = (u A (x, t), u A (x, t), ■ ■ ■ ,u A (x,t)) 
converging in (Ll oc (R} x IR + )) r to an entropy solution U e (x,t) = (u\,u\,--- ,u e r ) of the 
reaction-hyperbolic system fli.il) with initial data Uq(x). Furthermore, the solution fulfills 
the following estimates 

+ 0O 

(5.4) \U € {x,t)\ < \ U j ~ U j-i\ f° r almost aU M)> 

j=-oo 

+ 0O 



(5.5) TV(U*(;t))< \Uj-Ui-i\ 



(5.6) \\U*(.,t) - < C(hxp(- Xmm ^ tl} )\\Q°\\» +l)\t- t^. 
for all t, t\ > 0. 

In the framework of BV-solutions, the zero-relaxation limit can be very easily dis- 
cussed. In fact, the standard argument in [3] proves that the embedding of L°°(R + x 
R) n L°°(R+,BV(R)) n Lip(R + ,L 1 {R)) into Lj oc (R + x R) is compact. On the other 
hand, the estimates in (15. 4ft — ( 15. 6 j) show that {U e } 6> o lies in a bounded subset of 
L°°(R + xl)n L°°(R + , BV(R)) n Lip(R +1 L 1 (R)) by assuming 

(5.7) \\Q°\W = 0- 

Namely, the initial data are assumed to be in equilibrium. Thus, we have 
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Theorem 5.3. Under the conditions of Theorem \5.2\ and the equilibrium assumption 



( [5. ty , there exist a bounded measurable function U*(x, t) = {u\(x, t), u^x, t), ■ • ■ , u*(x, t)) 
and a subsequence(denoted in the same way) of set {U e (x,t) = {u\,u%, • ■ ■ , uf)} such that 
as e — > 0, 

U e (x,t) -»■ U*(x,t) in (LL( Rl x R+ )) r - 

Moreover, the function U*(x,t) is a weak entropy solution to the Cauchy problem (ji.gp 
with initial data Uq(x) and satisfies 

+oo 

\U*(x,t)\< \ U j~ U j-i\ for almost all (x,t), 

j=-oo 

+oo 

TV(U*(-,t))< \ U ?- U U\> 



||cr(-,t)-cr(.,ti)|| £ i <c|t-tn 

Remark 5.1. The fact that U*(x,t) satisfies the entropy conditions for the equilibrium 
system (jl.2p follows from Lemmas 14. II and 14. 21 Indeed, because S(U) in (14.61) is symmetric 
and non-positive, the term r]jj(U n+1 )Q(U n+1 ) in Lemma [4.21 is 

Vu{U n+1 )Q{U n+1 ) = nu{U n+l )S{U n+l )r]u{U n+l ) < 

if T) is chosen to be that constructed in Lemma 14.11 

Remark 5.2. Without the equilibrium assumption in (15. 7p . U e can only converge to U* 
for t > but not up to t = 0, because 

exp(-Amin{t,t 1 }/e) < — — , 6 

eAmmji, ii) 

in (15. 6p . Indeed, without the equilibrium assumption, initial boundary-layers occur. See 
also [9]. 
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